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ABSTRACT: Detrended, modelled first leaf dates for 856 sites across North America for the period 1900-2008 are used to 
examine how the El Nino Southern Oscillation (ENSO) and the Pacific Decadal Oscillation (PDO) separately and together 
might influence the timing of spring. Although spring (mean March through April) ENSO and PDO signals are apparent in 
first leaf dates, the signals are not statistically significant (at a 95% confidence level ( p <0.05)) for most sites. The most 
significant ENSO/PDO signal in first leaf dates occurs for El Nino and positive PDO conditions. An analysis of the spatial 
distributions of first leaf dates for separate and combined ENSO/PDO conditions features a northwest- southeast dipole 
that is significantly (at p <0.05) different than the distributions for neutral conditions. The nature of the teleconnection 
between Pacific SST’s and first leaf dates is evident in comparable composites for detrended sea level pressure (SLP) in 
the spring months. During positive ENSO/PDO, there is an anomalous flow of warm air from the southwestern US into 
the northwestern US and an anomalous northeasterly flow of cold air from polar regions into the eastern and southeastern 
US. These flow patterns are reversed during negative ENSO/PDO. Although the magnitudes of first leaf date departures 
are not necessarily significantly related to ENSO and PDO, the spatial patterns of departures are significantly related to 
ENSO and PDO. These significant relations and the long-lived persistence of SSTs provide a potential tool for forecasting 
the tendencies for first leaf dates to be early or late. Copyright © 2011 Royal Meteorological Society 
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1. Introduction 

The timing of phenological events (seasonal changes in 
plant and animal activity) are sensitive to climatic vari- 
ability and change (Running and Hunt, 1993; Myneni 
et al., 1997; Schwartz, 1997, 1998; Menzel and Fabian, 
1999; Sparks and Crick, 1999; Schwartz and Reiter, 2000; 
Cay an et al., 2001; Parmesan and Yohe, 2003; Menzel 
et al., 2006; Schwartz et al., 2006). Long-term, synop- 
tic phenological observations are rare in North Amer- 
ica, but certain phenological traits of select species can 
be modelled from climate data and then mapped at 
the continental scale. For example, the Spring Indices 
‘suite of metrics’ (Schwartz and Reiter, 2000) use mod- 
els of lilac ( Syringa chinensis ‘Red Rothomagensis’) 
and honeysuckle ( Lonicera tatarica ‘Arnold Red’, and 
L. korolkowii ‘Zabeli’), validated with continental- scale 


* Correspondence to: G. J. McCabe, US Geological Survey, Denver 
Lederal Center, MS 412, Denver, CO 80225. 

E-mail: gmccabe@usgs.gov 

^ The contributions of these authors to this article were prepared as part 
of their official duties as United States Lederal Government employees. 


networks of actual phenological observations, as proxies 
for the occurrence of first leaf, first bloom, and other 
temperature- sensitive responses to climate. Further, these 
Spring Indices (SI) provide an effective approach for 
processing climate data into metrics that are broadly rep- 
resentative of many deciduous plant species in temperate 
regions, and thus provide a long-term (1900-2008) and 
spatially widespread dataset of relevant phenoclimatic 
indicators (Schwartz et al., 2006). 

Analyses of time series of SI first leaf and first 
bloom dates for sites across North America have iden- 
tified regional differences in the variability and long- 
term trends of these phenological indicators and have 
indicated long-term trends to an earlier onset of spring 
associated with warming temperatures, especially for the 
northwestern US (Schwartz and Reiter, 2000). As might 
be expected, these trends are evident in monthly aver- 
ages of springtime temperatures and actual phenological 
observations made for lilac and honeysuckle since 1956 
(Cay an et al., 2001). The ultimate causes of interannual 
variations, secular trends, and their regional differences, 
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Figure 1. (A) location of modelled first leaf date sites used in the study, (B) number of sites with first leaf dates each year from 1900-2008, 
(C) 2° grid cells, and (D) mean first leaf date (day-of-the-year with 1 January being day 1) for 2° grid cells for the period 1900-2008. 


however, are not well understood, and their predictability 
in both the near and long term remains unexplored. 

Spatiotemporal variability in the timing of spring arises 
from complex climatic variations that are synthesized 
in indices of large-scale sea-surface temperature (SST), 
or atmospheric pressure patterns, in important centres 
of action and their regional teleconnections. Examples 
include El Nino-Southern Oscillation (ENSO), the North 
Atlantic Oscillation (NAO), the Pacific North American 
(PNA) pattern, the Arctic Oscillation (AO) or Northern 
Annular Mode (NAM), the Pacific Decadal Oscillation 
(PDO), and the Atlantic Multidecadal Oscillation (AMO). 
Because they can summarize complex and sequenced 
weather and climatic conditions into a single value, phe- 
nologists have started to employ these indices in explain- 
ing phenophase variations for both plants and animals 
(Cay an et al., 2001; Stenseth et al., 2003; Hallet et al., 
2004; Cook et al., 2005; Gordo, 2007; Ault et al., 2011). 
Also, seasonal persistence in these geophysical indices, 
or lags in their teleconnections, could inform long-lead 
forecasts of the onset of spring and related phenomena. 

The ENSO and the PDO have been identified as impor- 
tant modes of global climate variability that influence the 
variability of temperature and precipitation for various 
regions around the world (Kiladis and Diaz, 1989; Man- 
tua and Hare, 2002). There have been many studies of the 
influences of ENSO and PDO on North American climate 
(Cayan and Peterson, 1989; Redmond and Koch, 1991; 
Kahya and Dracup, 1993; Latif and Barnett, 1996; Man- 
tua et al., 1997; Gershunov and Barnett, 1998; McCabe 
and Dettinger, 1999; McCabe and Muller, 2002; McCabe 
and Wolock, 2008). In addition, several studies have 
investigated how ENSO and PDO interact and modify 
the expected ENSO hydroclimatic influences in the US 
(Gershunov and Barnett, 1998; McCabe and Dettinger, 
1999; Gutzler et al., 2002; Kurtzman and Scanlon, 2007). 

As ENSO and PDO are important modes of the global 
climate system and have substantial influences on the 
climate of North America, this study addresses the effects 


of ENSO and PDO on the onset of spring in North 
America as represented by modelled first leaf dates. The 
objectives of this study are to improve understanding of 
the climate factors that influence the variability of the 
onset of spring in North America, and to determine the 
possibility to forecast the onset of spring using ENSO and 
PDO indices (i.e. spring NIN03.4 SSTs and spring PDO). 


2. Methods and data 

A dataset of modelled first leaf and first bloom dates 
has been developed by Schwartz and Reiter (2000) and 
extended by Schwartz et al. (2006). This dataset includes 
first leaf and first bloom dates for 856 sites (Figure 1(A)) 
varying in number from 291 sites in 1900 to 829 sites in 
1967 (Figure 1(B)). These first leaf dates were computed 
using a SI model that is calibrated and validated using 
observations of first leaf and first bloom dates of lilac 
and honeysuckle. The first leaf dates used in this study 
are a model-derived indicator of the onset of spring that 
integrates a variety of springtime climatological variables 
(Schwartz and Reiter, 2000; Schwartz et al., 2006; Ault 
et al., 2011). Details on the development of the Sis are 
given in Schwartz et al. (2006) and Ault et al. (2011); 
also, Table I has a summary of the steps used to compute 
the SI first leaf dates. Among other things, our particular 
variant of the SI model considers chill requirements, the 
amount of winter chilling required by a plant to terminate 
dormancy and ensure adequate bud break the following 
spring. This becomes irrelevant in the warm subtropics, 
so we did not calculate SI for stations south of 34 °N 
(Figure 1(A)). 

The variability in the number of sites with data is a 
result of the availability of the necessary climate data 
(e.g. daily maximum and minimum temperature data) for 
each site. Comparison of the first leaf and first bloom 
dates (not shown) indicates substantial similarity between 
the two datasets. We chose first leaf for our analyses 
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Table I. Summary of how the Spring Indices first leaf date is calculated. 


Step 1: Calculate the day-of-the-year (DOY) dates (after 1 October) when 1375, 1250, and 1250 (for the clone lilac Syringa 
chinensis ‘Red Rothomagensis’, clone honeysuckle Lonicera tatarica ‘Arnold Red’ (HS ‘AR’), and clone honeysuckle L. 
korolkowii ‘Zabeli’ (HS ‘ZB’), sub-models respectively) chilling hours (base temperature 7.2 °C) are reached. Average the 
three DOY dates to produce Spring Indices Composite Chill Date. 

Step 2: Determine the DOY dates when the following equations are valid (for each sub-model): 

Lilac: 1000 = (3.306 x MDS0) + (SYNOP x 13.878) + (0.201 x DDE2) + (0.153 x DD57) 

HS ‘AR’: 1000 = (4.266 x MDS0) + (SYNOP x 20.899) + (0.248 x DD57) 

HS “ZB”: 1000 = (2.802 x MDS0) + (SYNOP x 21.433) + (0.266 x DDE2) 

Where: 

MDSO = number of days since sub-model chilling satisfaction date 

SYNOP = number of high-energy synoptic events that have occurred since sub-model chilling satisfaction date (three-day 
cumulative degree-hour accumulations higher than 637 degree-hours (base temp. = 0.6 °C)) 

DD57 = growing degree-hour accumulation (base temp. = 0.6 °C) over the 3 days that were 5-7 days prior to the current date 

DDE2 = growing degree-hour accumulation (base temp. = 0.6 °C] over the 3 days that were 0-2 days prior to the current date 

Step 3: Adjust the raw DOY dates determined in Step 2 to account for geographical and temporal bias in the number of 
high-energy synoptic events, using the following equations: 

Lilac: Adjusted date = 0.358 — 0.234 x RAW 

HS ‘AR’: Adjusted date = 1.579 - 0.225 x RAW 
HS ‘ZB’: Adjusted date = -0.204 - 0.109 x RAW 
Where: 

RAW = raw DOY dates for each sub-model 

Step 4: Average the DOY dates from the three sub-models to produce Spring Indices First Leaf Date 


because it may be most relevant to large-scale ecosystem 
processes (e.g. vegetation greenup), and it occurs earlier 
in the year when the teleconnections with ENSO/PDO 
are strongest. 

Currently, there is no one perfect index for the timing 
of spring that applies at the continental scale. The 
advantages of the SI are twofold: first, the SI is queued 
to the start of the growing season, and necessarily means 
the same thing year after year. That it is tied to lilac and 
honeysuckle is unimportant; these linkages just provide 
proof of concept. The SI is analogous to a drought 
index; a drought index may not be able to characterize 
a specific tree species’ response to moisture anomalies, 
but the drought index does mean something ‘absolute’ in 
terms of local climate variability through time and across 
space. The second advantage of the SI is that it only 
responds to the atmosphere; meaning that it is free of 
local and, thus, noisy biological and ecological effects. 
This is important because it maximizes the chances of 
identifying important relations between large-scale modes 
of climate variability and the index itself; an advantage 
other plant-based or data-driven indices (such as green- 
up) might not have. 

To provide century-long time series of first leaf dates, 
missing data for individual sites were filled using a 
regression approach. For each year, a regression was 
developed between first leaf dates (the dependent vari- 
able) for sites with data and the latitude and longitude 
for these sites (the independent variables). The regression 
was then used to estimate first leaf dates for sites with 
missing data. The distribution of coefficients of determi- 
nation for the regressions for each year indicated reliable 
regression models. The coefficients of determination had 
a median value of 0.70, and a 25th percentile of 0.64 and 
a 75th percentile of 0.74. 


Because the density of modelled first leaf data is higher 
in eastern North America than in western North America, 
one additional step was performed to provide a more even 
spatial distribution of first leaf dates and to increase the 
smoothness of first leaf date spatial patterns. This step 
involved gridding the first leaf data for the 856 sites to 
2° x 2° grid cells (2° grid cells) (Figure 1(C)). A 2° x 2° 
grid resolution provides a large enough number of grid 
cells for robust statistical analysis and allows most of 
the grid cells to include multiple sites of first leaf data. 
Grid cells were retained if at least one site was located 
within a grid cell. This resulted in a set of 2° grid cells 
numbering 242. Separate principal components analyses 
(not shown) of (1) the point data and (2) the 2° gridded 
data provided nearly the exact same results. This result 
indicates that little or no information regarding first leaf 
date variability was lost by gridding the data to the 2° 
grids. 

Figure 1(D) illustrates the mean first leaf date for the 
period 1900-2008. Mean first leaf dates primarily occur 
during the months of March through May. Because the 
first leaf dates for most sites occur during the months 
of March through May, we focused on these months as 
the period of the analyses (Figure 1(D)). Additionally, 
the time series of first leaf dates for each 2° grid cell 
were converted to departures from the respective long- 
term (1900-2008) mean first leaf dates for each grid 
cell to make it easier to identify positive and negative 
departures through time and across space. 

Trends in first leaf dates indicate earlier first leaf 
dates for a large part of the western US and southwest- 
ern Canada (Figure 2). These trends to earlier first leaf 
dates are consistent with increases in spring temperatures 
in these regions (Cayan et al., 2001; Lu et al., 2005). 
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Figure 2. Trends in modelled first leaf dates during the 1900-2008 
period. The black squares indicate trends that are statistically significant 
at p < 0.05. 

Because trends can influence analyses of temporal vari- 
ability in data time series, the linear trends in the 2° 
gridded time series of first leaf dates were removed. 

NIN03.4 SSTs are a useful index of ENSO (Tren- 
berth, 1997). For this study, monthly NIN03.4 SSTs 
were computed from the Kaplan monthly SST dataset 
(http://www.esrl.noaa.gov/psd/data/gridded/data.kaplan_ 
sst.html). Monthly SSTs averaged for the region from 5 °S 
to 5 °N and from 170 °W to 120 °W were used to represent 
NIN03.4 SSTs. A time series of monthly PDO values 
was obtained from http://jisao.washington.edu/pdo. The 
monthly values of NIN03.4 SSTs and PDO were aver- 
aged for the months of March through May (MAM) to 
generate times series of MAM NIN03.4 SSTs and MAM 
PDO values for the years 1900-2008. Similar to the first 
leaf date time series, linear trends in MAM NIN03.4 and 
MAM PDO time series were removed for the analyses in 
this study. 

The analyses in this study primarily involve computing 
(and mapping) mean departures (in days) of detrended 
mean first leaf dates for combinations of positive and 
negative detrended MAM NIN03.4 and detrended MAM 
PDO conditions. Positive and negative MAM NIN03.4 
and MAM PDO conditions were determined based on 
terciles computed for the detrended MAM NIN03.4 and 
detrended MAM PDO time series (1900-2008). Thus, 
the upper (lower) terciles were used to identify positive 
(negative) MAM NIN03.4 and MAM PDO conditions. 
Neutral MAM NIN03.4 and MAM PDO conditions 
were identified as years when both the NIN03.4 and 
PDO values were in the middle tercile of the respective 
distributions. Terciles were used instead of quartiles to 
allow enough years to be included in each category for 
statistical analysis. 

For each site and sets of years for combinations of 
positive and negative MAM NIN03.4 and MAM PDO 
conditions, mean departures (in days) of first leaf dates 
were computed. The mean departures were subsequently 
mapped and the maps were examined and compared to 
identify influences of ENSO and PDO on the variability 
of first leaf dates. The number of observations (years) 
used for neutral March through May (MAM) NIN03.4 


Table II. Number of observations (years) used for neutral 
MAM NIN03.4 and PDO conditions and separate and com- 
bined positive (+) and negative (— ) MAM NIN03.4 and MAM 
PDO conditions. 


NIN03.4 and/or PDO Condition 

Number of Observations 

neutral 

16 

+NIN03.4 

36 

— NIN03.4 

35 

+PDO 

36 

-PDO 

35 

+NIN03.4, +PDO 

20 

— NIN03.4, +PDO 

5 

+NIN03.4, -PDO 

5 

— NIN03.4, -PDO 

17 


and PDO conditions and separate and combined positive 
and negative MAM NIN03.4 and MAM PDO conditions 
are listed in Table II. 

Student’s r-tests were used to identify mean departures 
that are significantly different (at a 95% confidence level, 
p <0.05) from neutral conditions. In addition, statisti- 
cal tests were performed to determine if the number of 
statistically significant departures was significantly differ- 
ent from what would be expected (field significance test) 
(Eivezey and Chen, 1983; McCabe and Dettinger, 1999). 
The determination of whether the number of negative 
and positive departures was statistically significant (at 
p <0.05) was accomplished through a Monte Carlo anal- 
ysis in which 1000 ‘pseudo-NIN03.4’ and ‘pseudo-PDO’ 
time series were generated using a stochastic model. The 
simulated time series were 109 years in length (the same 
length as the record of first leaf data), and the simulated 
time series preserved the lag-1 autocorrelations in the 
NIN03.4 and PDO time series as well as the underlying 
correlation between NIN03.4 and PDO (Newman et al ., 
2003). 

The 1000 simulated time series of detrended MAM 
NIN03.4 and PDO were then used to select sets of 
years to compute mean departures of 2° gridded first leaf 
dates for each combination of positive/negative MAM 
NIN03.4 and MAM PDO conditions. For each of the 
1000 simulations the numbers of negative and positive 
significant departures (at p <0.05) were computed. From 
these values, the 95th percentile values of the number 
of significant negative and positive departures were 
computed and used as thresholds of statistical significance 
(field significance) for the numbers of significant negative 
and positive departures resulting from the original first 
leaf data for separate and combined MAM NIN03.4 and 
MAM PDO conditions. 


3. Results and discussion 

Figure 3 illustrates mean departures of first leaf dates for 
positive and negative MAM NIN03.4 SSTs and positive 
and negative MAM PDO. For these composites, NIN03.4 
and PDO effects are considered separately. Given that 
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Figure 3. Mean departure of modelled first leaf date for positive (+) and negative (— ) MAM NIN03.4 SSTs and MAM PDO during 1900-2008. 
The black squares indicate departures that are statistically significant at p <0.05. 


MAM NIN03.4 and MAM PDO are related SST indices, 
the spatial patterns of mean departures of first leaf dates 
are expectedly similar for both SST indices. For positive 
MAM NIN03.4 (Figure 3(A)) and positive MAM PDO 
(Figure 3(B)), departures of first leaf dates are negative 
(earlier first leaf dates) in the northwestern US and south- 
western Canada and are positive (later first leaf dates) 
in the southeastern US. For negative MAM NIN03.4 
(Figure 3(C)) and negative MAM PDO (Figure 3(D)), 
departures of first leaf dates are positive in the north- 
western US and southwestern Canada and negative in 
the southeastern US. These patterns of departures of first 
leaf dates are consistent with patterns of temperature 
anomalies in North America associated with NIN03.4 
and PDO (Kiladis and Diaz, 1989; Redmond and Koch, 
1991; Mantua and Hare, 2002). When NIN03.4 and PDO 
are positive (negative), positive (negative) temperature 
anomalies occur in the northwestern US and, conse- 
quently, negative (positive) temperature anomalies occur 
in the southeastern US. 

For positive MAM PDO (Figure 3(B)) the mean depar- 
tures of first leaf dates are significantly ( p <0.05) dif- 
ferent from neutral conditions for a number of grid 
cells, whereas there are only a few grid cells with 
significant departures for positive and negative MAM 
NIN03.4 (Figures 3(A) and (C)), and negative MAM 
PDO (Figure 3(D)). An analysis of the number of grid 
cells with statistically significant (at p <0.05) differ- 
ences from neutral conditions indicates that the number 
of significant differences for none of the composite maps 
exceed the expected value (determined using a field sig- 
nificance test and a significance level of p <0.05). 

Gershunov and Barnett (1998) suggest that when 
NIN03.4 SSTs and PDO are in similar phases that the 
effects of these SST patterns on North American climate 
are amplified (i.e. constructive phases). To examine 
these effects, mean departures of first leaf dates for 


combinations of positive/negative MAM NIN03.4 and 
positive/negative MAM PDO were computed. 

For combined positive MAM NIN03.4 and positive 
MAM PDO conditions, the mean departures of first leaf 
dates are negative (earlier first leaf dates) in the north- 
western US and southwestern Canada (Figure 4(A)). Pos- 
itive departures (later first leaf dates) primarily occur in 
the southeastern US. This spatial pattern of first leaf date 
departures is consistent with the effects of El Nino (and 
positive PDO) conditions on temperatures across North 
America (Redmond and Koch, 1991; Mantua and Hare; 
2002). During El Nino (and positive PDO) temperatures 
in the northwestern US and southwestern Canada are gen- 
erally warmer than average (resulting in earlier first leaf 
dates), and temperatures in the southeastern US are gen- 
erally cooler than average (resulting in later first leaf 
dates) (Redmond and Koch, 1991; Mantua and Hare, 
2002 ). 

In contrast, during combined negative MAM NIN03.4 
and negative MAM PDO conditions, the spatial pattern 
of first leaf date departures indicates positive departures 
in the northwestern US (and southwestern Canada) and 
negative in the southeastern US (Figure 4(D)). This 
pattern of departures of first leaf dates is negatively 
correlated with the pattern for combined positive MAM 
NIN03.4 and positive MAM PDO (Pearson correlation 
coefficient (r) = —0.32, p <0.01). 

For combined positive MAM NIN03.4 and negative 
MAM PDO conditions (Figure 4(C)), the spatial pattern 
of first leaf departures indicates a mixture of departures; 
whereas, for combined negative MAM NIN03.4 and pos- 
itive MAM PDO (Figure 4(B)), the first leaf departures 
are negative for most of the western US and southwestern 
Canada, and positive for most of the eastern US. 

The pattern of first leaf departures for combined 
positive MAM NIN03.4 and positive MAM PDO 
(Figure 4(A)) appears to be the strongest of the patterns 
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Figure 4. Mean departure of modelled first leaf date (in days) for combinations of positive (+) and negative (— ) MAM NIN03.4 SSTs and 
MAM PDO during 1900-2008. The black squares indicate departures that are statistically significant at p <0.05. 


of combined positive MAM NIN03.4 and negative MAM 
PDO conditions, with 20 grid cells indicating significant 
negative departures and 7 grid cells with significant pos- 
itive departures. However, these numbers of significant 
positive and negative departures is less than the expected 
number for the number of grid cells, using a 95% confi- 
dence level. 

3.1. Differences across space 

Although the mean departures of first leaf dates did 
not indicate many statistically significant departures (at 
p <0.05) for separate and combined MAM NIN03.4 
and MAM PDO conditions, the effects of NIN03.4 and 
PDO could be seen in the patterns of departures and 
these patterns are consistent with the effects of NIN03.4 
and PDO on North American temperature (Redmond and 
Koch, 1991; Mantua and Hare, 2002). The effects of 
NIN03.4 and PDO on first leaf date departures also can 
be seen in the frequency distributions of first leaf date 
departures for separate and combined MAM NIN03.4 
and MAM PDO conditions. Figure 5 illustrates frequency 
distributions of first leaf date departures (expressed as 
probability density functions (pdfs)) for separate and 
combined MAM NIN03.4 and MAM PDO conditions). 
The pdfs of first leaf date departures for NIN03.4 and 
PDO conditions are flatter and wider than is the pdf 
for neutral conditions. Flatter and wider pdfs indicate 
fewer values near the centre of the range of values, 
and more values toward the tails of the distribution. 
The flatter and wider pdfs for NIN03.4 and PDO 
conditions, compared with the pdf for neutral conditions, 
suggests an effect of NIN03.4 and PDO on first leaf date 
departures. Additionally, the pdfs for positive NIN03.4 
and positive PDO (Figure 5(A)) are flatter and wider 
than are the pdfs for negative NIN03.4 and negative 
PDO. For combinations of NIN03.4 and PDO the pdfs 
for positive NIN03.4 and positive PDO, and negative 


NIN03.4 and positive PDO are flatter and wider than 
any of the pdfs for other combined effects (Figure 5(B)). 
These results suggest that positive PDO conditions have 
a stronger effect on first leaf departures than does 
NIN03.4. 

The Index of Agreement ( d ) was used to compare the 
distributions of first leaf date departures for separate and 
combined MAM NIN03.4 and MAM PDO conditions 
with the distribution of first leaf date departures for 
neutral conditions (Willmott, 1981, 1984; Legates and 
McCabe, 1999). The Index of Agreement is often used 
to compare two time series (e.g. measured and modelled 
time series), but it also can be used to compare samples 
of paired data. For our application, d is computed as 

^ (Sh-Shf 

d=l. N ^ ,=1 _ = — (1) 

T. A\S2i - Sl\ - \Sh - Sl\) 2 

*-^l = 1 

where N is the number of sites, Sli is the sample 1 
value for site i, S2 t is the sample 2 value for site i, 
and SI is the mean for all sites in sample 1. The Index 
of Agreement ranges from 0 to 1, with a value of 1 
indicating perfect agreement between two samples. The 
Index of Agreement represents an improvement over the 
coefficient of determination (r 2 ) for sample comparisons 
because it is sensitive to differences in both means and 
variances of two samples (Legates and McCabe, 1999). 

The d statistics for all of the comparisons of the 
first leaf date distributions for separate and com- 
bined NIN03.4 and PDO conditions with the distri- 
bution for neutral NIN03.4 and PDO are near zero 
(Table III). These results indicate that the first leaf 
date distributions for separate and combined NIN03.4 
and PDO conditions are substantially different from 
the distribution for neutral conditions. These results 
indicate that the effects of NIN03.4 and PDO are 
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Figure 5. Probability density functions of the distributions of mean 
departures of modelled first leaf dates for (A) positive (+) and 
negative (-) MAM NIN03.4 SSTs (NIN03.4) and MAM PDO, and 
(B) combinations of positive and negative MAM NIN03.4 and MAM 
PDO. The probability density function of mean first leaf date departures 
for neutral MAM NIN03.4 and MAM PDO conditions is plotted (black 
line) as a reference for the other distributions. 


Table III. Index of Agreement ( d ) for comparisons of dis- 
tributions of mean first leaf date departures for separate and 
combined positive (+) and negative (— ) MAM NIN03.4 and 
MAM PDO conditions with the distribution of mean first leaf 
date departures for neutral MAM NIN03.4 and MAM PDO 
conditions. 


d 


+NIN03.4 

0.03 

— NIN03.4 

0.09 

+PDO 

0.06 

-PDO 

0.10 

+NIN03.4, +PDO 

0.08 

— NIN03.4, +PDO 

0.09 

+NIN03.4, -PDO 

0.13 

— NIN03.4, -PDO 

0.14 


strongly manifested in the spatial patterns of first leaf 
dates. 

3.2. Relations with atmospheric circulation 

An examination of sea level pressure (SLP) anomalies 
for the years included in each of the composites for 


positive (negative) MAM NIN03.4 and positive (nega- 
tive) MAM PDO help explain the patterns of the first 
leaf departures (Figure 6). As with other datasets, the 
SLP anomalies were detrended before analyses were per- 
formed. Since the patterns of mean first leaf departures 
for the constructive combinations of MAM NIN03.4 and 
MAM PDO conditions (i.e. positive NIN03.4 and pos- 
itive PDO, and negative NIN03.4 and negative PDO) 
included the largest number of observations and indi- 
cated the strongest effects of NIN03.4 and PDO on first 
leaf dates, SLP anomalies for these two conditions are 
examined (Figure 6). 

For combined positive MAM NIN03.4 and posi- 
tive MAM PDO conditions (Figure 6(A)), positive SLP 
anomalies are centred over central Canada and the north- 
central US; whereas, negative SLP anomalies occur over 
the eastern North Pacific Ocean and extend into west- 
ern North America (i.e. a strengthened Aleutian Low). 
Negative SLP anomalies also occur over the western 
North Atlantic Ocean and extend into the southeastern 
US. These SLP anomalies suggest an anomalous flow of 
warm air from the southwestern US into the northwest- 
ern US and an anomalous northeasterly flow of cold air 
from polar regions into the eastern and southeastern US. 
These SLP anomalies and the associated anomalous flows 
of air are consistent with warmer-than-average spring 
temperatures (early first leaf dates) in the northwestern 
US and southwestern Canada, and cooler- than- average 
spring temperatures (late first leaf dates) in the south- 
eastern US. 



-IS -ft -G -3 Q 3 G 9 IS 

Mean Departure of First Leaf Date (days) 


Figure 6. Mean departure of modelled first leaf date (in days) for 

(A) positive MAM NIN03.4 SSTs and positive MAM PDO, and 

(B) negative MAM NIN03.4 SSTs and negative MAM PDO, during 
1900-2008, and associated SLP anomalies (in hectoPascals). The solid 
lines indicate positive SLP anomalies and dashed lines indicate negative 

SLP anomalies. The contour interval is 0.2 hectoPascals. 
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In contrast, for combined negative MAM NIN03.4 and 
negative MAM PDO conditions (Figure 6(B)) negative 
SLP anomalies occur over Canada and the north-central 
US and positive SLP anomalies occur over the North 
Pacific and North Atlantic Oceans. This pattern of SLP 
anomalies is the opposite of the pattern of SLP anomalies 
for positive NIN03.4 and positive PDO conditions (r = 
—0.68, p <0.01). These SLP anomalies indicate an 
anomalous northwesterly flow of cool air (and late 
first leaf departures) into southwestern Canada and the 
northwestern US, and an anomalous southwesterly flow 
of warm air (resulting in negative first leaf departures) 
into the southeastern US. 

The first leaf date departures in western North America 
and the associated SLP anomalies for positive NIN03.4 
and PDO (Figure 6(A)) conditions and negative NIN03.4 
and PDO conditions (Figure 6(B)) are consistent with 
results for ENSO and PDO relations with ice duration 
on lakes and rivers in Canada given by Bonsai et al ., 
2006, who found that ENSO and PDO had substantial 
effects on ice duration. They found that positive phases 
of ENSO and PDO are associated with shorter ice 
durations consisting of later freeze-ups and earlier break- 
ups. In contrast, negative phases of ENSO and PDO 
were associated with longer ice duration over much 
of Canada. These results are consistent with previous 
analyses that showed that positive ENSO and PDO 
conditions are associated with warmer temperatures over 
much of Canada, particularly in western Canada during 
winter and spring (e.g. Bonsai et al ., 2001). Bonsai et al. 
(2006) suggest that the higher temperatures over western 
Canada and resultant shorter ice durations during positive 
ENSO and PDO conditions are mainly attributable to 
a strengthened Aleutian Low and associated warm air 
advection into much of Canada (Figure 6(A)). 

The analysis of SLP anomalies associated with the spa- 
tial patterns of first leaf date departures indicates that 
ENSO and PDO variability influence the timing of first 
leaf in North America through effects of SLP patterns 
on temperature.These analyses provide a physical expla- 
nation for the variability of departures of first leaf dates 
for combinations of NIN03.4 and PDO conditions. In 
addition, although the magnitudes of first leaf departures 
are not strongly related to NIN03.4 and PDO variability, 
the spatial patterns of first leaf dates and the tendencies 
to be early or late are substantially related to NIN03.4 
and PDO. 


3.3. Potential for forecasts of first leaf date departures 

The relations between first leaf dates and NIN03.4 and 
PDO, and the persistence of SSTs (Table IV) suggest 
potential for forecasts of the tendency of first leaf dates to 
be early or late. The relations of NIN03.4 and PDO with 
SLP anomalies (Figure 6) also suggest the possibility of 
using SLPs to forecast the tendency for early or late first 
leaf dates. However, the atmosphere does not have as 


Table IV. Correlations of mean three-month average NIN03.4 
SSTs, PDO, and the PNA index with mean average MAM 
values. 



NIN03.4 

PDO 

PNA 

OND 

0.75 

0.57 

0.01 

NDJ 

0.77 

0.64 

0.03 

DJF 

0.80 

0.76 

0.03 

JFM 

0.84 

0.87 

0.59 

FMA 

0.87 

0.96 

0.81 

MAM 

1.00 

1.00 

1.00 


much persistence as do SSTs. For example, the PNA 
index characterizes large-scale atmospheric circulation 
over North America (Wallace and Gutzler, 1981; Yarnal 
and Diaz, 1986; Yarnal and Leathers, 1988) and has been 
recognized as a primary mode of atmospheric circulation 
variability over North America (Barnston and Livezey, 
1987; Leathers and Palecki, 1992). The PNA serves as 
a measure of the relative amplitude of the stationary 
waves over North America and has proven to be useful 
to understand climate variability over on time scales 
ranging from months to decades (Leathers and Palecki, 
1992). Additionally, the PNA is related to variability in 
NIN03.4 and PDO (Yu and Zwiers, 2007). 

Correlations between MAM PNA and the PNA aver- 
aged for three-month periods before the spring sea- 
son show that the persistence of PNA is short-lived 
(Table IV), whereas the persistence of NIN03.4 and PDO 
is large and long-lived. Thus, SSTs likely provide the best 
opportunity for potential forecasts of first leaf dates. 

To explore the potential for forecasting first leaf 
date departures, time series of average NIN03.4 and 
PDO for the OND period were used to compute the 
risk of early first leaf dates. Risks were computed for 
the constructive combinations of NIN03.4 and PDO 
(i.e. positive NIN03.4 and positive PDO conditions, 
and negative NIN03.4 and negative PDO conditions). 
There are 18 years during 1900-2008 with positive 
OND NIN03.4 and positive OND PDO conditions, and 
19 years with negative OND NIN03.4 and negative 
OND PDO conditions. Risk of early first leaf dates was 
computed as the number of years indicating first leaf dates 
in the lowest percentile during 1900-2008 at each 2° grid 
cell (expressed as a percentage of the total number of 
years for the constructive combinations of OND NIN03.4 
and OND PDO). These results suggest that there may be 
potential to forecast at least the spatial pattern of first 
leaf date departures when NIN03.4 SSTs and PDO are 
in positive phases. 

Results indicate that for positive OND NIN03.4 and 
positive OND PDO conditions the risk of early first leaf 
dates is greater than 50% for much of the northwestern 
US and southwestern Canada (Figure 7(A)). For these 
NIN03.4 and PDO conditions, the risk of early first leaf 
dates is lowest (less than 10%) in the southeastern US. 
This pattern of risks of early first leaf dates suggests that 
useful forecasts of early first leaf dates may be possible 
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Figure 7. Risk of early modelled first leaf date for (A) positive OND 
NIN03.4 and positive OND PDO conditions, and (B) negative OND 
NIN03.4 and negative OND PDO conditions during 1900-2008. 


for several months ahead of the spring season when 
positive NIN03.4 and positive PDO conditions exist. 

For negative OND NIN03.4 and negative OND PDO 
conditions the pattern of risks of first leaf dates indicates 
low risks for most locations across the northern US and 
southern Canada, and higher risks across the southern 
US (Figure 7(B)). The highest risks of early first leaf 
dates for these conditions are not as high as the highest 
risks indicated for positive OND NIN03.4 and positive 
OND PDO conditions (Figure 7(A)). This result is not 
surprising given that the persistence of negative NIN03.4 
and negative PDO conditions does not appear to be 
as strong as is the persistence of positive NIN03.4 
and positive PDO conditions. Thus, the forecasts using 
negative NIN03.4 and negative PDO conditions likely 
are less reliable than those made using positive NIN03.4 
and positive PDO conditions. 


4. Conclusions 

Time series of SI first leaf dates, an index of the onset 
of spring, for 856 sites across North America for the 
period 1900-2008 were analysed to examine the effects 
of ENSO and PDO variability on the timing of spring. 
Results indicate noticeable (although not necessarily 
statistically significant) temporal effects of ENSO and 
PDO on the timing of first leaf dates in the northwestern 
US and southwestern Canada. During El Nino and 
positive PDO conditions, first leaf dates are generally 
earlier in the northwestern US and later in the eastern US. 
With La Nina and negative PDO conditions, the opposite 
pattern of first leaf date departures generally occurs. 


Although the temporal effects of ENSO and PDO on 
first leaf dates is not necessarily statistically significant, 
analyses of differences in the spatial patterns of first 
leaf dates associated with ENSO and PDO indicate 
highly statistically significant differences from neutral 
conditions for all separate and combined conditions of 
NIN03.4 and PDO. In addition, the spatial patterns of 
first leaf date departures can be explained by differences 
in SLP anomalies driven by NIN03.4 and PDO. Results 
also suggest that the persistence of NIN03.4 and PDO 
when in positive phases may be useful to forecast spatial 
patterns of first leaf date departures. 
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